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Experimental  Study  of  a Tunnel 's 
Collapse  Criterion 


by 

Leif  N.  Persen 


Introduction. 

This  presentation  is  an  attempt  to  gather  the  informa- 
tion obtained  through  experimental  investigations  of  what  in  general  terms 
might  be  called  a cavity's  "collapse  criterion"  as  far  as  dynamic  loading 
is  concerned.  The  experiments  to  be  considered  are  the  ones  sponsored  by 
the  Bundesministerium  der  Verteidigung,  West  Germany  and  carried  out  by 
A/ S NORCONSULT,  Oslo,  Norway  in  the  period  1968-74. 

It  should  be  emphasized  right  from  the  beginning,  that  the 
collapse  criterion  in  the  sense  used  here  is  not  a single-valued  physical 
quantity  even  for  a given  tunnel  in  a specific  rock.  A tunnel  may  have 
several  "collapse  criteriae"  according  to  the  definition  of  what  should 
be  considered  "a  collapse".  The  purpose  of  the  tunnel  will  here  be  one 
of  the  many  factors  to  be  considered.  A tunnel  may  serve  such  a purpose 
that  it  will  continue  to  function  until  a more  or  less  complete  cave-in 
occurs.  On  the  other  hand  its  purpose  may  be  to  protect  highly  sensitive 
objects  from  shocks  caused  by  oncoming  stress  waves,  in  which  case  its 
protective  function  may  cease  long  before  the  tunnel  as  such  is  endangered. 
The  collapse  criterion  is  thus  here  closely  linked  to  the  cessation  of 
the  tunnels  functioning  according  to  its  purpose. 

The  very  wide  definition  given  above  and  the  almost  infi- 
nite variety  of  possible  situations  which  may  occur  made  it  necessary  to 
limit  the  scope  of  the  investigation  to  a few  situations  of  practical 
importance.  First  a distinction  was  made  between  the  lined  and  the  unlined 
("raw")  tunnel.  Then  only  a certain  class  of  linings  were  picked  to  be 
tested.  Finally  limitations  were  put  upon  the  ratio  between  the  tunnel 
diameter  and  the  length  of  the  oncoming  stress  wave.  Within  this  rather 
restricted  framework  only  limited  information  on  the  interaction  between 
the  tunnel  and  the  stress  wave  was  sought. 
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1.  The  "raw"  tunnel. 


The  ability  of  a "raw"  tunnel  to  withstand  the  influence 
of  an  oncoming  shock  wave  without  being  damaged,  depends  on  several  condi- 
tions which  may  be  of  a rather  vague  character.  One  such  condition  will 
be  the  way  in  which  the  tunnel  was  originally  excavated.  As  indicated  in 


Fig.l  Dislocation  zones  around  two  similar  tunnels  under 
different  conditions . 

Fig.l  the  dimensions  of  the  zones  of  dislocated  rock  material  around  the 
tunnel  may  vary  considerably  depending  upon  the  method  of  blasting  used 
to  excavate  it  as  well  as  upon  the  original  condition  of  the  rock  before 
excavation. 

Another  such  condition  will  be  the  extent  to  which  static 
stresses  are  present  in  the  rock  material  adjacent  to  the  tunnel.  These 
static  stresses  may  change  with  time,  and  even  though  these  changes  may 
be  small  and  very  slow,  the  stage  may  be  reached  at  which  the  tunnel  may 
collapse  without  any  additional  dynamic  loading. 

These  two  conditions  show  how  difficult  it  will  be  to 
establish  valid  "collapse  criteriae"  for  "raw"  tunnels.  They  also  illu- 
strate the  difficulties  involved  in  a theoretical  approach  to  the  problem 
because  such  qualities  as  "degree  of  dislocation"  is  not  easily  quanti- 
sized  as  part  of  a description  of  the  situation  around  the  tunnel. 

The  experimental  approach  to  the  problem  consisted  of  an 
attempt  to  perform  an  experiment  which,  when  executed  on  location,  might 
give  an  idea  of  the  collapse  criterion  of  a given  tunnel.  The  idea  was 
to  excavate  a test  tunnel  using  the  same  procedure  as  that  used  when  the 
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real  tunnel  was  excavated.  Because  one  operates  on  location  one  might  hope 
to  have  the  same  static  stresses  as  those  present  in  the  real  tunnel 
neglecting  the  long  term  rearrangement  which  the  stress  at  the  real  tunnel 
may  have  undergone.  Because  of  the  way  in  which  the  test  tunnel  is  exca- 
vated one  might  also  hope  to  have  dislocation  zones  which  are  similar 
in  the  two  cases.  The  test  tunnel  will  preferably  be  made  smaller  than 
the  real  tunnel,  and  consequently  the  question  of  a model  law  will  arise. 

The  idea  behind  the  tests  was  to  create  an  inward  mov- 
ing concentric  shock  wave  around  a circular  cylindrical  tunnel.  One  could 
thereby  hope  to  take  advantage  of  the  focusing  effect  of  the  shock  wave 
[l]  so  that  one  could  create  large  enough  stresses  for  impending  "collapse" 
at  the  test  tunnel  without  endangering  the  real  tunnels  in  the  neighbor- 
hood. This  idea  led  to  an  experimental  set-up  illustrated  in  Figs. 2 and  3. 

A circular  tunnel  was  excavated  with  a vertical  axis.  The  radius  of  the 
tunnel  was  R3  = 0.70  [m]  . On  a concentric  circle  of  radius  Ri,  = 1 . 70  [m] 
holes  for  pick-ups  based  on  strain  gauges  were  drilled.  By  means  of  these 
pick-ups  the  inward  moving  stress  wave  can  be  monitored.  Holes  for  the 
charges  were  drilled  on  circles  with  radiae  Ri  = 3.10  [m]  and  R2  = 2.U0  [m] . 
Cylindrical  charges  were  used  as  illustrated  in  Fig. 3 . On  the  walls  of  the 
tunnel  were  mounted  accelerometers  in  the  region  where  the  shock  wave 
was  expected  to  be  "plane".  Fig. 2 shows  the  real  location  of  the  charges, 
the  pick-ups  and  the  accelerometers  in  the  experiment.  A similar  experi- 


Fig.  2 Geometry  of  the  experiment  (top  view) 


ment  was  also  made  to  see  if  the  focusing  effect  of  the  shock  wave  was 
achieved.  In  that  case  the  geometry  shown  in  Fig. 2 was  kept,  the  tunnel 
was  however  not  excavated  but  the  accelerometers  were  replaced  by  ordi- 
nary strain  gauge  pick-ups  placed  in  holes  at  the  same  location. 
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Fig.  3 Geometry  of  the  experiment  (side  view) 


The  result  of  these  tests  can  be  summed  up  as  follows: 

1.  As  could  be  expected,  the  stress  waves  from  each  source  did  not 
merge  into  one  single  nice  concentric  shock  wave.  Analysis  of  the 
the  shock  wave  showed  that  it  would  always  have  to  be  considered 
as  the  sum  of  the  cylindrical  charges  superimposed  on  each  other. 
This  was  however  sufficient  for  the  purpose. 

2.  The  focusing  effect  of  the  geometry  was  counteracted  by  the 
attenuation  in  such  a way  that  the  peak  value  of  the  shock  wave 
remained  constant  as  it  moved  inward.  This  was  considered  satis- 
factory. 

3.  The  collapse  criterion  was  to  be  determined  as  that 
value  of  the  oncoming  shock  wave  for  which  spalling  of  the 
walls  started  to  occur.  The  charges  were  steadily  increased 
and  it  turned  out  that  no  definite  value  could  be  found  below 
which  spalling  did  not  occur.  A certain  amount  of  spalling 
could  be  traced  even  at  very  low  levels.  Thus  it  was  decided 
to  use  the  mass  of  the  spall  as  a measure  for  "collapse",  and 

a limiting  value  of  360  [atm]  for  the  permissible  peak  value  of 
the  oncoming  shock  wave  was  finally  established  as  the  "collapse 
criterion"  for  this  particular  tunnel  in  that  particular  rock. 
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One  would  have  to  take  account  of  the  fact,  that  by  steadily 
increasing  the  charges,  the  repeated  detonations  may  alter  the 
rock's  properties  as  a wave  transmitting  medium. 


The  walls  of  the  tunnel  were  painted  in  a color  contrasting  that 
of  the  rock  so  that  the  places  where  spalling  had  occurred  could 
be  easily  found.  In  this  way  the  mounted  picture  of  the  wall 
after  the  event  shown  in  Fig. 4 was  obtained. 


Fig. 4 A mount  of  single  photos  of  the  tunnel  walls  after  spalling 
areas  indicate  places  where  spalling  has  occurred. 


The  results  obtained  in  these  experiments  needed  additional 
support,  and  supplementing  experiments  were  made  in  an  entirely  different 
type  of  rock.  One  of  the  contentions  to  be  tested  was  among  others  the  assumed 
geat  influence  on  the  criterion  of  the  method  used  to  excavate  the  tunnel. 

Two  of  the  test  tunnels  were  therefore  drilled,  leaving  the  neighboring 
rock  as  good  as  untouched,  whereas  the  other  two  were  excavated  by  cautious 
blasting.  The  geometry  of  the  experiment  is  illustrated  in  Figs. 5, 6 and  7. 

The  test  tunnels  I,  II,  III  and  IV  extend  from  the  four  walls  of  the  chamber 
A in  Fig. 5 . The  tunnels  have  a radius  R 


indicated  in  Fig. 6 . 
On  concentric  circles  were  drilled  holes  for  the  pick-ups  (A1  - A8) , and 
for  the  charges  (B1  - B16) . On  the  tunnel  walls  were  mounted  accelerometers 
(Cl  - C6) . The  length  of  the  charges  were  4m  as  shown  in  Fig. 7 where  also 
the  relative  postion  of  the  pick-ups  and  the  accelerometers  is  shown. 

The  charges  were  increased  in  steps  until  spalling  took 
place  to  such  an  extent  that  the  tunnel  was  considered  useless  as  shelter. 
The  value  of  the  maximum  amplitude  of  the  shock  wave  for  which  this  occurred 
was  taken  as  the  collapse  criterion.  It  is  however  clear  that  this  value 
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Fig.  7 Geometry  of  the  experiment.  ( Cut  along 
the  axis  of  the  tunnel.  ) 

may  be  influenced  by  the  repeated  detonations.  Therefore  the  charges 
were  raised  from  30  to  120  [g/m]  by  tunnel  II  whereas  it  was  decreased 
from  120  to  30  [g/m]  by  tunnel  IV.  The  effect  on  the  mass  of  the  spall 
was  striking.  The  two  tunnels  could  be  assumed  to  be  rather  equal  and 
thus  could  be  expected  to  behave  similarly.  The  mass  of  the  spall  by 
120  [g/m]  was  far  less  when  this  was  the  first  detonation  than  when 
one  had  gradually  increased  the  charge  to  the  same  value. 

The  influence  of  the  method  used  to  excavate  the  tunnels 
was  also  clearly  brought  out.  The  mass  of  the  spalls  in  tunnels  I and  III 
which  were  drilled,  was  much  less  than  in  the  other  two  which  were  blasted. 
If  it  were  not  for  two  discontinuity  surfaces  intersecting  the  drilled 
tunnels,  where  the  only  spalling  took  place,  spalling  would  most  probably 
not  have  occurred  at  all  in  these  tunnels  at  the  same  level  at  which  the 
others  had  to  be  considered  "collapsed". 

The  pick-ups  (based  on  strain  gauges)  and  the  accelero- 
meters were  used  to  monitor  the  shock  wave.  The  measurements  on  the  tunnel 
walls  corresponded  well  with  those  made  in  the  rock. 
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The  results  of  these  experiments  beyond  what  has  already 
been  mentioned  are  perhaps  best  illustrated  in  the  following  diagrams. 
From  the  signals  obtained  by  the  accelerometers  at  the  walls  of  the 
tunnels,  the  radial  velocities  at  the  walls  could  be  obtained  as  func- 
tions of  time.  Fixing  the  attention  on  the  peak  values  one  found  that 
the  same  charges  measured  in  grams  per  meter  ({/)  did  not  always  create 
tlie  same  wall  velocities  in  the  tunnels.  The  peak  values  would  also 
vary  from  one  accelerometer  to  another.  One  may  however  take  the  mean 
peak  value  of  the  wall  velocities  and  correlate  with  the  charge  mag- 
nitude as  has  been  done  in  Fig. 8 . One  will  then  find  that  a rather 
nice  correlation  existed  in  the  case  of  tunnels  II,  III  and  IV. 


0.5  1.0  1.5  2.0 


Mean  peak  velocity  [m/  s] 

Fig.  8 Correlation  between  charge  magnitude  W and  the  mean  peak  velo- 
city created  at  the  tunnel  walls  of  the  different  tunnels. 

It  is  however  easily  recognized  that  tunnel  I deviates  from  the  rest 
and  thus  the  mean  value  of  the  peak  radial  velocities  measured  by 
means  of  the  accelerometers  at  various  locations  for  one  shot  should 
be  correlated  with  the  magnitude  of  the  spalling  caused  by  the  shot. 
This  is  done  in  Fig. 9 . The  magnitude  of  the  spalling  is  measured  as 
the  weight  of  the  spalls  (A/)  and  Fig. 9 reveals  that  now  tunnel  III 
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Mean  peak  velocity  Qn/s] 


Fig. 9 Mass  of  spall  M shown  as  function  of  the  mean  peak  velocity 
measured  at  the  wall  for  each  shot 


seem  to  deviate  from  the  others.  It  should  however  be  kept  in  mind 
that  the  tunnels  I and  III  were  drilled,  and  that  spalling  in  these 
tunnels  took  place  only  at  the  locations  were  these  tunnels  were  inter- 
sected by  the  discontinuity  surface  in  the  rock.  The  result  in  Fig. 9 
for  these  tunnels  are  thus  not  representative  for  a tunnel  drilled 
in  the  undisturbed  rock. 

The  limiting  value  of  M below  which  the  spall  is  considered 
to  be  of  tolerable  magnitude  has  to  be  judged,  and  will  sometimes  de- 
pend on  other  considerations  than  the  functioning  of  the  tunnel  as 
[ ’ such.  In  Fig. 9 this  limit  is  chosen  such  that  the  limiting  value  is 

^ , reached  when  the  peak  wall  velocity  is  1.1*  [m/s]  . This  corresponds  to 

a maximum  amplitude  of  the  oncoming  shock  wave  of  103  [atm] , which  is 
less  than  1/3  of  the  value  obtained  in  the  previous  case. 

One  may  now  sum  up  the  results  of  these  investigations 
in  the  following  points  valid  for  raw  tunnels: 

k 1.  It  has  been  verified  that  the  method  used  to  excavate  the  tunnel 

will  greatly  influence  the  collapse  criterion  of  the  tunnel. 

2.  This  means  however  in  view  of  the  long  terms  changes  in  the 
state  of  static  stress  around  the  tunnel,  that  also  the  age 
of  the  tunnel  will  influence  the  collapse  criterion  of  the  tunnel. 
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3.  Using  test  tunnels  on  location  and  performing  the  described 
type  of  tests  on  them  may  for  a given  location  give  an  indica- 
tion of  the  collapse  criterion  in  the  sense  used  here.  However, 
the  uncertainties  of  such  a procedure  are  so  great,  that  one  must 
contemplate  very  carefully  if  the  expenditures  are  warranted. 
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2.  The  1 ined  tunnel . 

The  experiment- ai  investigations  of  the  collapse  criterion 
of  lined  tunnels  to  be  reported  on  here  were  all  carried  out  in  the 
same  rock.  Because  one  had  to  limit  the  number  of  possible  types  of 
linings  to  a managable  number,  only  the  types  exhibited  in  Fig. 10  were 
considered.  These  represent  the  normal  types  of  linings  used  when  more 
sophisticated  arrangements  are  not  needed. 


Fig. 10  Lined  tunnels  to  be  tested,  a)  Concrete  lining  b)  Reinforced 
concrete  lining  a)  Reinforced  concrete  lining  with  anchors 


Because  the  tunnels  must  be  expected  to  be  subjected  to 
the  shock  waves  from  sherically  symmetric  sources,  the  geometry  of  the 
experimental  set-up  was  chosen  such  that  the  test  tunnels  would  be  sub- 
jected to  spherical  shock  waves.  It  was  also  decided,  that  because  one 
would  in  most  cases  want  to  test  tunnels  in  scales  less  than  1:1,  the 
arrangement  should  preferably  be  such  as  to  give  information  on  the 
influence  of  the  size  of  the  tunnel  radius.  What  could  be  called  a 
model  law  for  such  experiments  was  envisaged. 

In  addition  to  the  questions  outlined  above,  main  empha- 
sis was  to  be  placed  upon  an  evaluation  of  the  importance  of  anchors, 
i.e.  how  much  would  the  carrying  capacity  of  the  tunnel  be  increased 
by  the  introduction  of  anchors  as  compared  with  the  same  lining  with- 
out anchors.  A similar  consideration  was  also  to  be  made  to  determine 
the  importance  of  the  reinforcement. 

These  considerations  led  to  a geometry  of  the  experiment 
which  is  sketched  in  Fig. 11.  In  the  floor  of  a large  chamber  were 
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Fig. 11  Sketch  of  the  geometry  of  the  experimental  set-up. 


excavated  seven  vertical  tunnels  numbered  1 through  7.  Tunnels  1,  4,  5 
and  7 had  radiae  of  R = 1.7  [it]  , tunnels  2,  3 and  6 had  radiae  of 
R = 0.85  [m]  . Tunnels  1 and  2 were  unlined  ("raw"),  tunnels  3 and  4 
were  lined  as  b),  tunnels  6 and  7 were  lined  as  a)  and  tunnel  5 as  c) , 

a),  b)  and  c)  referring  to  the  type  of  lining  shown  in  Fig.  10. 

These  tunnels  were  to  be  subjected  to  shock  waves  created 
in  I and  II,  two  wells  into  which  the  charges  were  to  be  detonated. 

One  could  not,  for  practical  reasons,  make  these  wells  so  small  that 
the  charges  to  be  used  could  be  conceived  of  as  completely  confined. 

The  charges  were  therefore  detonated  hanging  free  in  the  center  of 
the  wells  which  were  filled  with  water.  In  this  way  it  was  also  hoped 
to  create  identical  shock  waves  in  I and  II  using  the  water  to  transmit 
the  input  pulse  to  the  rock.  For  testing  the  measuring  devises  small 
confined  charges  were  detonated  in  III,  IV  and  V,  which  were  small  holes 
drilled  in  the  rock.  The  distance  from  the  main  charge  in  I or  II  to  the 

tunnels  was  5 [m]  as  shown  in  Fig. 11  . However,  the  unlined  tunnels  1 

and  2 were  placed  at  a distance  of  12  [m]  from  II,  because  their  capacity 
to  withstand  the  effect  of  the  shock  wave  was  estimated  to  be  far  less 
than  the  lined  tunnels'. 

The  tunnels  were  10  [m]  deep  as  shown  in  Fig. 12  . They 
were  lined  to  a height  of  6 [m]  above  the  floor.  The  charges  were  always 
placed  at  a depth  of  7 [m]  below  the  top  surface,  at  which  depth  the 
measuring  sensors  were  placed.  These  sensors  consisted  of  accelerome- 
ters attached  to  the  tunnel  walls  as  shown  in  Fig. 13.  In  addition  a row 
of  strain  gauge  based  pick-ups  were  placed  between  the  two  wells  for 
the  charges.  The  signals  from  these  were  used  to  get  information  on 


> vi 

1 1 

1 

1 

Fig. 13  Top  view  of  the  tunnels  showing  the  positions  of  the  accelerometers. 
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the  shock  wave  created  in  each  case  by  the  explosion.  From  the  accelero- 
meter signals  information  could  be  obtained  on  the  motion  of  the  walls 
and  the  maximum  amplitude  of  the  shock  wave  when  reaching  the  tunnel  could 
be  determined.  , 

The  intention  was  to  use  charges  large  enough  to  bring 
the  tunnel  wall  to  failure.  The  maximum  amplitude  of  the  stress  wave 
as  it  reached  the  tunnel  could  then  be  used  as  a collapse  criterion. 

One  would  however  have  to  try  to  get  failure  at  the  first  attempt  be- 
cause the  previous  experiments  had  shown  a marked  influence  of  repeated 
explosions.  One  started  out  with  1 and  2 kg  charges  in  III,  IV  and  V 
for  calibration  purposes.  It  was  however  decided  also  to  detonate  small- 
er charges  in  I and  II  for  the  same  purpose,  but  because  of  the  way  in 
which  the  shock  wave  was  created  the  charges  were  here  10  kg. 

The  first  real  test  was  the  detonation  of  50  kg  TNT  in  II. 
This  was  done  to  fix  at  least  one  lower  level  which  was  not  dangerous 
to  the  tunnels.  No  visible  damage  was  done  to  any  tunnel  on  this  occa- 
sion. Then  200  kg  TNT  was  detonated  in  I and  thereafter  in  II.  In  these 
cases  heavy  damage  was  inflicted  on  all  tunnels  except  the  unlined 
tunnels  1 and  2,  which  remained  intact  due  to  the  longer  distance  from 
the  source. 

Some  conclusions  were  drawn  from  inspection  of  the  damages: 

1.  Comparison  between  tunnels  of  different  sizes  but  with  the  same 
type  of  lining  (tunnels  3 and  4 and  to  some  extent  also  tunnels 
6 and  7)  showed  that  the  yield  mechanisms  were  basicly  the  same 
in  both  cases.  The  damages  seemed  however  to  be  greater  in  the 
lager  tunnels,  an  observation  which  will  be  given  some  thoughts 
later . 

2.  Comparison  between  tunnels  of  the  same  size  but  with  different 
types  of  lining  (tunnels  4,  5 and  7)  showed  that  the  damages 
inflicted  on  tunnel  7 (concrete  lining  without  reinforcement) 
were  much  larger  than  for  the  other  two.  Tunnel  7 collapsed 
completely  as  is  seen  from  the  photo  in  Fig. 14  . It  was  con- 
cluded that  in  most  cases  would  this  kind  of  lining  contribute 
very  little  to  increasing  the  tunnel's  capasity  to  withstand 
the  influence  of  the  shock  wave. 

3.  The  yield  mecl  anism  by  lined  tunnels  under  these  conditions  may 

be  contemplated  by  observing  the  fracture  of  the  walls  of  tunnel  4 
shown  in  Fig. 15. 

4.  The  fracture  of  tunnel  4 is  similar  to  that  of  tunnel  5.  The  differ- 
ence in  the  lining  of  the  two  tunnels  consists  of  the  anchors 

and  one  may  conclude  that  although  the  anchors  may  have  contri- 
buted to  change  the  stress  pattern  slightly,  the  collapse  cri- 
terion is  for  all  practical  purposes  the  same.  The  influence  of 
the  anchors  is  therefore  marginal. 
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Fig. 14  View  of  tunnel  7 after  the  200  kg  shot.  The 
tunnel  has  collapsed  completely  and  only  a 
small  portion  of  the  original  wall  can  still 
be  identified. 


Fig. IS  View  of  tunnel  4 after  the  200  kg  shot.  T 
yield  occurred  at  two  distict  gene  rati ces 
The  yield  mechanism  is  shown  skematical'y 
in  Fig.  7 6 


TUNNEL  No.  4 


Fig. 18  Mount  of  photos  of  the  wall  of  tunnel  4 after  the  200  kg  shot 
The  yield  occurred  as  in  tunnel  5 but  the  yield  zone  is  here 
located  at  0 - 54  . The  generatices  shown  indicate  the  loca- 

tion of  the  Accelerometers. 


The  yield  mechanism  suggested  by  the  fracture  zones  of 

tunnels  4 and  5 is  given  as  a shear  fracture  in  Fig. 16  . This  occurs 

along  generatrices  as  clearly  indicated  by  the  photo  mounts  in  Figs. 

17  and  18.  One  observes  of  course  spalling  also  elsewhere  on  the 

walls  whereby  the  reinforcement  is  uncovered,  but  the  real  fracture 

of  the  lining  occurs  at  an  angle  8 . This  angle  is  different  in  the 
o o . ^ 

two  tunnels  (75  in  4 and  54  in  5)  and  this  may  be  attributed  to  the 
influence  of  the  anchors.  In  a later  section  this  will  be  dicussed 
in  context  with  an  attempt  to  explain  the  fracture  in  detail. 

The  smaller  tunnels  3 and  6 were  also  brought  to  collapse 
The  damage  inflicted  upon  them  did  not  seem  to  be  as  great  as  for  the 
same  type  of  tunnel  with  a larger  radius.  Fig. 19  shows  tunnel  6 after 


Fig. 19  View  of  tunnel  6 after  the  200  kg  shot.  The 
fracture  takes  place  along  generatices  as 
in  the  case  of  tunnel  4 

the  200  kg  shot.  The  tunnel  should  be  compared  with  tunnel  7,  a tunnel 
with  the  same  lining  but  with  a larger  radius.  Tunnel  7 collapsed  com- 
pletely, tunnel  6 showed  collapse  with  damage  comparable  to  tunnel  4 
and  with  fracture  occurring  along  generatrices.  It  should  be  noted, 
that  the  thickness  of  the  lining  was  the  same  in  tunnels  6 and  7,  and 
that  the  thicknes/radius  ratios  were  different  in  the  two  cases.  This 
will  later  be  jsed  to  explain  the  difference  in  the  damage. 

Fig. 20  shows  tunnel  3 after  the  200  kg  shot.  This  tunnel 
had  the  same  lining  as  tunnel  4,  it  collapsed  in  the  same  way  as  this 
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Fig.  20  View  of  tunnel  3 after  the  200  kg  shot.  The 
fracture  takes  place  along  generatrices  as 
in  the  case  of  tunnel  4 


tunnel  and  again  the  damage  was  observed  to  be  smaller  for  the  tunnel 
with  the  smaller  radius.  All  comments  to  tunnel  6 apply  also  to  this 
case . 

This  concludes  the  remarks  on  the  damages  as  such  and 
represent  also  the  answer  to  some  of  the  questions  originally  posed. 
However,  the  instrumentation  gave  additional  information  of  quanti- 
tative nature  which  could  be  used  to  supplement  the  information  ob- 
tained so  far.  Such  information  would  concern  the  wall  velocities 
created  by  the  shock  wave  etc.  For  that  purpose  the  undisturbed 
tunnels  1 and  2 were  equipped  with  a lining  of  the  type  b)  in  Fig. 10 
and  subjected  to  the  influence  of  shock  waves  created  in  a well  as 
described  here.  For  further  supplementation  the  lining  was  doubled 
in  thickness  to  see  if  the  effect  of  the  thickness  of  the  lining 
could  be  detected  in  the  results.  Again  the  tunnels  were  subjected 
to  the  influence  of  shock  waves,  whereby  in  both  cases  the  maximum 
amplitude  of  the  shock  waves  was  kept  low  enough  to  avoid  damage  to 
the  linings.  The  results  of  these  experiments  are  described  in  the 
following  section. 


J 


20 


3.  The  wal 1 motion. 


The  experiments  described  so  far  had  the  aim  to  provide 
answers  to  rather  limited  questions,  but  as  mentioned  one  found  that 
answers  of  greater  generality  might  perhaps  be  found.  Hence  the  ex- 
tention  of  the  experiments  was  approved. 

One  finds  that  cases  may  occur  in  which  the  contents 
hidden  in  the  shelter  and  their  sensitivity  determine  the  extent 
to  which  the  tunnel  may  be  subjected  to  the  influence  of  shock  waves. 
Thus  cases  may  occur  where  the  "collapse  criterion"  is  an  upper  bound 
on  the  maximum  permissible  amplitude  of  the  shock  wave  at  the  tunnel 
which  is  far  below  the  value  for  which  the  tunnel  as  such  is  endangered. 
The  way  in  which  the  contents  will  experience  the  shock  wave  will  de- 
pend on  the  arrangement  within  the  tunnel,  but  in  all  cases  the  motion 
of  the  tunnel  walls  during  the  period  when  the  tunnel  is  engulfed  by 
the  shock  wave  will  be  of  significant  importance. 


WELL  FOR  CHARGES 


ACCELEROMETERS 


Fig.  21  Geometry  of  the  supplementary  experiments  showing  the 
two  test  tunnels  1 and  2 , the  well  for  the  charges  and 
the  row  of  strain  gauge  based  pick-ups  for  monitoring  the 
stress  wave. 
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The  geometry  of  the  supplementary  experiments  is  shown 
in  Fig. 21  . It  is  noticed  that  the  tunnels  are  now  equipped  with  accel- 
erometers mainly  on  that  side  which  faces  the  oncoming  shock  wave.  The 
reason  for  that  is  that  in  the  first  experiments  very  little  move- 
ment of  the  "back  walls"  was  recorded  even  in  the  cases  where  the 
"front  walls"  were  heavily  damaged.  It  should  be  noticed  that  this 
indicates  a special  type  of  collapse  , and  that  thus  the  results  are 
limited  to  this  situation. 

Before  going  into  detail  on  the  results  the  following  con- 
siderations are  useful.  The  interaction  between  a shock  wave  and  a tun- 
nel will  to  some  extent  be  characterized  by  the  ratio  D/L  between  the 
diameter  D of  the  tunnel  and  the  length  L of  the  oncoming  shock  wave. 

When  D » L the  interaction  will  be  similar  to  the  reflection  from  a 
plane  wall,  and  the  maximum  amplitude  of  the  stress  wave  may  have  been 
considerably  attenuated  when  travelling  a distance  equal  to  the  dia- 
meter of  the  tunnel.  Usually  this  situation  is  not  of  any  great  prac- 
tical importance.  When  D « L the  variation  in  the  stress  across  the 
tunnel  is  not  significant  compared  with  the  stress  level  in  the  wave. 

One  has  a quasi-static  case  with  a stress  distribution  around  the  tun- 
nel which  may  be  judged  from  static  considerations.  This  may  be  the 
case  when  the  stress  wave  originates  from  a nuclear  source.  When  D - L 
the  situation  occurs  which  the  experiments  are  supposed  to  cover,  and 
which  is  rather  difficult  to  handle  theoretically.  In  the  tests  carried 
out  in  the  experiments  the  ratio  varied  within  the  range  1.36  < L/D  < 5.67. 

The  scope  of  the  investigation  was  to  establish  relations 
which  would  permit  the  prediction  of  the  motion  of  the  tunnel  walls  when 
the  oncoming  shock  wave  was  specified.  Such  a specification  can  be  given 
by  3 quantities  if  the  oncoming  wave  is  considered  to  be  mainly  a tri- 
angular wave  as  shown  in  Fig. 22  . The  maximum  amplitude  A is  measured  in 
[m/s]  whereas  the  wave  length  L and  the  rising  length  A L are  measured 


Fig.  22  Simplified  form:  of  the  oncoming  shock  wave 


in  [ s j and  converted  to  lengths  by  means  of  the  signal  velocity  c • of 
the  rock.  This  really  means  that  the  shock  wave  is  characterized  by  its 
velocity  field. 


One  should  bear  in  mind  that  measurements  of  wave  lengths 
are  very  uncertain  and  that  great  care  must  be  exercised  when  deciding 
whether  or  not  irrelevant  influences  are  present  in  the  signals.  As  an 
example  Fig. 23  shows  the  signals  for  shots  4-8  with  the  triangular 
wave  drawn  in.  The  "tail"  of  the  original  signal  is  influenced  by  re- 
flections and  has  therefore  been  neglected  in  deciding  the  triangular 
replacement  wave.  f^/^. 


SHOT  4 


SHOT  5 


SHOT  6 


SHOT  7 


SHOT  8 


23.  Wave  forms  by  shots  4-8  with  their  triangular  replan  m 
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With  the  three  input  data  A,  L and  A L in  addition  to 

the  signal  velocity  c of  the  rock  the  motion  of  the  tunnel  walls 

3 

can  be  predicted  in  the  following  sense: 

1.  Just  as  attention  is  fixed  only  at  the  first  positive  phase 

of  the  oncoming  shock  wave,  only  the  first  positive  phase 

of  the  radial  wall  velocity  V is  considered. 

v 

2.  The  first  positive  phase  of  V is  conceived  of  as  having  a 
triangular  form. 

3.  The  wall  motion  will  then  be  given  by  four  quantities: 

V given  as  a function  of  the  location  8,  0 being 

p 1TICLX  ° * ° 

J the  angle  at  the  center  of  the  tunnel  between 

the  direction  to  the  source  and  the  direction 

to  the  location. 

t given  as  the  time  at  which  motion  at  the  location 

a n 

0 starts. 

t given  as  the  rising  time,  i.e.  the  time  needed 

for  V to  reach  its  maximum  value  at  location  0. 
r 

A t given  as  the  duration  of  the  positive  phase  of 

the  radial  wall  velocity  at  location  0 . 

The  aim  is  to  relate  these  four  quantities  to  the  three  parameters 

of  the  shock  wave. 

The  arrival  time  t : 

a 


Fig. 24  Geometry  of  a spherical 
shock  wave  engulfing  a 
cylindrical  cavity. 

The  arrival  time  t is  counted  from  the  time  when  the  shock 
wave  first  meets  the  cavity  until  it  has  reached  the  position  9 as  show’n 
in  Fig. 24.  From  purely  geometrical  considerations  one  obtains: 
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t R , , / r /?,  t ' 1 

— - #{V  1 + 2 T + (T)2  (1  - cose;  " 2 } 
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This  may  be  replaced  by  the  much  simpler  expression  when  R,  » R : 

K 


- (1  - cosQ  ) 


The  maximum  amplitude  V 

r,  max 

The  maximum  amplitude  V of  the  radial  component 

v,max 

of  the  wall  velocity  may  according  to  the  findings  of  the  investiga- 
tions be  expressed  as  follows: 


V = V exp(-(c tijpr)2)  (3.3) 

where  a is  a constant  expressing  the  influence  of  the  size  of  the 
tunnel  and  where  0 is  to  be  introduced  in  degrees  (°) . The  size  of 
the  tunnel  is  related  to  the  decay  of  the  stress  wave  as  it  travels 
a distance  comparable  to  the  radius  of  the  tunnel.  Let  be  the  non- 
dimensional  peak  stress  in  the  shock  wave  as  its  front  impinges  on 
the  tunnel.  Let  S ^ be  the  corresponding  value  when  the  wave  has 
travelled  so  far  that  its  front  has  reached  the  position  0 - 0.  Then 
a is  given  by: 


a = 1.  784  S /Snn  + 0.256 
o 90 


It  is  clear  that  in  this  way  not  only  the  size  of  the  shock  wave  as 
compared  with  the  tunnel,  but  also  the  "properties"  of  the  rock  as 
a wave  transmitting  medium  has  been  attempted  accounted  for. 

The  duration  time  At  ; 

The  duration  time  At  of  the  first  positive  phase  of  the 
radial  wall  velocity  may  be  expressed  as  follows: 


At  - A toexp{-((^)2) 


where  8 according  to  the  resuLts  seem  to  be  independant  of  the  tunnel 
size  and  consequently  equal  to  a constant: 

6 - 2.25  (3.6) 

The  quantity  A t is  the  duration  time  recorded  at  the  location  0 - 0} 
and  it  may  be  related  to  the  oncoming  shock  wave  as  follows: 

At  = — tanh(5R/L)  (3.7) 

s 

where  R is  equal  to  the  tunnel  radius  as  before  and  L is  the  wave 
length  of  the  shock  wave  in  the  sense  used  previously. 

The  rising  time  t : 

The  time  needed  for  the  velocity  at  a given  location  to 
increase  from  zero  to  its  first  peak  is  called  the  rising  time  t and 
may  be  given  as  a fraction  of  the  duration  time: 

t yA  t = 0.059  + 1.510(Q/180)2  (3.8) 

The  numerical  constants  in  this  expression  are  rather  uncertain.  They 
have  been  determined  from  information  from  the  first  experiments  in  the 
geometry  shown  in  Fig. 11,  and  Fig. 25  gives  an  impression  of  the  accuracy 


Fig.  26  The  ratio  t /At  plotted  as  function  of  0.  Results  from  experi- 
ments in  the  geometry  of  Fig.  21 


tJAt  = 0.  1477  + l.S60(Q/180)2 
t /At  - 0.1079  + 1. 5S2(Q/180)2 
t yAt  - 0.  1318  + 2.  078(Q/ 180)  2 
t /At  - 0.059  + 1.510(Q/180)2 


(Tunnel  2) 

(Tunnel  1) 

(Tunnel  2) 
(Tunnels  1 and  2) 


The  later  experiments  in  the  geometry  of  Fig. 21  gave  results  exhi- 
bited in  Fig. 26.  The  uncertainty  of  the  constants  is  clearly  demon- 
strated. It  should  be  noticed  that  this  uncertainty  is  immaterial 
for  the  determination  of  the  maximum  radial  displacement  of  the 
wall,  it  is  however  crucial  if  one  wants  to  determine  the  maximum 
radial  acceleration  of  the  wall. 


Fig. 27  Comparison  between  measured  wall  velocities  at  tunnt  l 7 I 

dotted  line)  and  the  predicted  velocities  based  on  the  equa- 
tions (3.1)  to  (3.8)  ( full  drawn  line) 


The  way  in  which  the  empirical  formulae  (3.1)  to  (3.8)  may 
be  used  is  exhibited  in  Figs. 27  and  28.  In  these  diagrams  the  recorded 
wall  velocities  are  shown  as  functions  of  time.  These  are  compared  with 
the  signals  which  can  be  predicted  based  on  the  given  formulae  provided 
the  oncoming  shock  wave  is  specified.  Fig.  27  shows  the  comparison  for 
the  50  kg  shot  against  tunnel  7 where  the  shock  wave  specification  has 
been  gathered  from  the  recordings  made  of  the  wave  propagation.  It 
should  be  noted  that  this  is  a tunnel  with  a large  radius,  and  with 
a lining  that  did  not  stand  up  to  the  200  kg  shot.  Fig. 28  shows  the 
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same  type  of  comparison  for  tunnel  3 (Fig. 11)  which  is  a tunnel  with 
only  half  the  radius  of  tunnel  7 and  with  a different  lining.  It  is 
seen  that  in  both  cases  the  predictions  occur  with  a satisfactory 
degree  of  accuracy  bearing  in  mind  the  great  uncertainties  of  such 
measurements.  It  seems  that  the  type  of  lining,  at  least  within  the 
variations  considered  here,  does  not  influence  the  wall  motion  as 
long  as  the  stresses  are  kept  on  a level  below  collapse. 

4.  A semi-static  approach. 

The  investigations  described  here  are  only  aimed 
at  gaining  an  understanding  of  the  collapse  of  a tunnel  when  this 
takes  place  as  a puncture  of  the  lining  during  the  passage  of  the 
first  phase  of  the  shock  wave.  Only  this  case  has  been  studied 
experimentally,  and  the  empirical  relations  obtained  are  limited  to 
such  cases.  However,  it  may  be  of  interest  to  attempt  a more  de- 
tailed study  of  how  collapse  is  brought  about  under  such  circum- 
s tances . 

If  one  compares  the  carrying  capacity  of  a "raw" 
(unlined)  tunnel  with  that  of  a lined  one,  it  seems  natural  to 
attribute  the  increased  capacity  of  the  latter  to  the  fact  that 
the  lining  is  carrying  the  burden  exclusively.  The  lining  will 
act  as  a shell  or  an  arch  which  during  the  passage  of  the  shock 
wave  carries  a rapidly  changing  load.  One  may  attempt  computing 
the  stresses  in  the  lining  by  a static  approach,  assuming  that^ 
at  each  position  of  the  shock  wave  relative  to  the  tunnel,  the 
lining  acts  as  an  arch  carrying  a static  load  which  however 
would  have  to  vary  from  one  position  to  another.  The  problem 
will  then  be  to  determine  the  load  to  be  applied. 

In  the  following  the  situation  will  be  examined 
when  the  front  of  the  shock  wave  passes  the  front  portion  of 
the  tunnel.  The  tunnel  is  supposed  to  have  a circular  cross  section, 
the  front  of  the  shock  wave  has  a given  slope  and  its  position 
during  the  passage  is  given  by  the  angle  9 as  shown  in  Fig. 29  . 

Two  positions  of  the  shock  wave  are  shown  in  Fig.  29  , and  the 
load  function  p giving  the  distribution  of  the  load  on  the  arch 
over  its  horizontal  projection  is  supposed  to  be  proportional  to 
the  amplitude  of  the  shock  wave  in  an  undisturbed  field.  This 
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Fig. 29  . The  front  of  the  shook  wave  at  two  positions 
0 = 0q  with  the  corresponding  load  function  p. 

then  leads  to  the  following  expression  for  the  load  function 


p - p (cosQ  - cosQ  ) 
r r o o 


(4.1) 


It  should  be  stressed  that  this  is  to  be  regarded  as  a guess,  and  that 
a similar  distribution 


P = Pr 


1 - cosQ 


sinQ 


(sinQ  - sinQ) 
o 


(4.2) 


II 


would  give  a straight  line  distribution  with  the  same  value  at  0 = 0 . 

o 

This  would  be  just  as  acceptable  a guess  for  the  load  function,  and 
consequently  also  this  possible  distribution  will  be  considered. 

It  is  further  assumed  that  the  lining  acts  as  an  arch 

which  is  built  in  at  its  ends  where  0 = ±0  . The  situation  will  then 

o 

be  as  sketched  in  Fig. 30  where  the  loading  situation  on  the  arch  is 
given.  H and  V are  the  horizontal  and  the  vertical  components  re- 
spectively of  the  force  transmitted  through  the  arch  at  its  built- 

in  end  and  M is  the  bending  moment  at  the  same  location.  At  an  ar- 

o 

bitrary  location  0 the  forces  transmitted  through  a cross  section 
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Fig.  30  . Sketch  of  the  arch  with  its  load  for  a given  position  of  the 
shock  wave.  The  positive  directions  of  the  forces  at  an  ar- 
bitrary cross  section  are  also  shown. 

of  the  arch  are  given  by  the  shearing  force  Q , the  axial  force  N and 
the  bending  moment  M.  These  quantities  have  been  calculated  as  func- 
tions of  6 for  each  position  9=9  of  the  shock  wave,  and  the  re- 

o 

suit  is  given  in  the  tables  in  the  Appendix.  The  load  function  used 
is  given  in  (4.1)  but  for  comparison  also  (4.2)  has  been  used.  The 
result  can  be  summed  up  through  Fig. 31  where  two  characteristic 
situations  are  exhibited.  The  distribution  of  the  axial  force  N, 
the  shearing  force  Q and  the  bending  moment  M over  the  arch  for  two 

different  positions,  9 = 54°  and  9 = 81°,  of  the  shock  wave  are 

o o 

shown.  The  load  funtion  is  given  in  (4.1)  but  in  the  latter  case 
also  the  results  for  the  load  function  in  (4.2)  are  shown  by  dotted 
lines.  It  is  noticed  that  the  difference  between  the  two  load  func- 
tions is  marginal  for  all  quantities  except  the  axial  force  V.  This 
quantity  exhibits  a maximum  value  which  does  not  occur  at  9 = 9^  as 
Q and  M do.  Comparing  these  results  to  the  experimentally  obtained 
yield  of  the  lining  (Figs. 17  and  18)  one  finds  that  the  position  at 
which  N exhibits  ist  maximum  value  corresponds  very  well  to  the  po- 
sitions at  which  yield  occurred.  Thus  one  may  take  this  as  an  encourag- 
ing indication  that  the  proposed  semi-static  approach  may  render 
sensible  results. 

Regarding  the  yield  mechanism  for  a reinforced  concrete 
structure  subjected  to  the  simultaneous  action  of  axial  - and  shear 
stress  as  well  as  bending  stresses  no  conclusive  investigation  seems 
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Fig.  31.  Axial  force  N}  shear  force  Q and  bending 
moment  M in  the  arch  (lining)  for  two 
different  positions  of  the  shock  wave. 

to  be  available  in  the  literature.  A recent  investigation  origi- 
nating from  the  activities  in  the  North  Sea  has  been  carried  out 
at  the  Norwegian  Institute  of  Technology  (NTH)  by  LENSCHOW  and 
HOFS0Y  [2].  The  complexity  of  the  situation  is  perhaps  best  brought 
out  by  a brief  account  of  their  results.  Fig. 32  is  a reproduction 
of  Fig. 5. 4 in  [2].  Is  shows  how  the  carrying  capacity  of  a rein- 
forced concrete  structure  depends  on  the  magnitude  of  axial-  and 
shear  forces  as  well  as  of  bending  moments,  and  indicates  the  type 
of  failure  under  different  circumstances  according  to  the  Norwegian 
Code  NS  3473.  It  is  however  clear  that  the  code  given  this  way  is 
valid  only  for  a given  relation  between  the  shear  force  and  the 
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LOAD  CAPASITY 

According  to  the  Norwegian  Code, 
NS  3473 


Fig  5.4 


N0=  frbh 
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X iNo 
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SHF.  A It  CAPASITY. 
(Compression  failure)  y I 


I MOMENT  CAPASITY^  I 


''SHEAR  CAPASITY 
(Diagonal  tension  failure) 


0.2  U / 
I / 


1 23456789 

Fig. 32  . Reproduction  of  Fig. 5. 4 of  [ 2] - 


bending  moment.  (The  notations  used  in  Fig. 32  and  later  in  Fig. 33 

are  explained  in  the  figure.  It  should  be  added  that  / and  f are 

J o J s 

the  design  strength  of  the  concrete  and  the  steel  respectively.  V 
is  the  shear  stress  corresponding  to  Q in  the  present  notations.) 

Fig. 33  is  a reproduction  of  Fig. 9. 8 of  [2]  which  gives 
a review  of  the  situation.  Not  only  is  the  Norwegian  Code  3473  shown, 
but  also  the  American  Code  ACI  and  the  European  Code  CEB  together 
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Analytical  model  compared  with  codes 
LOAD  CAPASnV 


1,20  0.60  OAO  0,30  0,24  0,20  0,17  0,13  0,13 


•-Tests; 


A.  ■ Aj*0,67_b!i 
s ‘ 100 


Vh'™ 

%■ 


0,82 


★ . MODEL  TEST 

(Related  to  the  higher 
concrete  strength) 

☆ > MODEL  TEST 

(Not  related  to  the  higher 
concrete  strength) 


Fig. 33  . Repvoduation  of  Fig. 9. 8 of  [2] 


with  their  own  results  and  their  analytical  model.  It  seems  clear 
that  these  results  indicate  a shear  failure  at  the  locations  ob- 
served if  the  results  from  the  semi-static  calculations  are  used. 
Thus,  one  is  again  encouraged  to  try  this  approach  for  the  deter- 
mination of  the  capacity  of  a lined  tunnel  to  withstand  the  influence 
of  a shock  wave.  Again  it  is  stressed  that  the  approach  only  is  aimed 
at  cases  where  the  failure  takes  place  as  a "puncture"  of  the  lining 
during  the  passage  of  the  first  phase  of  the  shock  wave. 


'i-  «%•* 


35 


5.  Acknowledgement 


The  present  investigation  has  been  sponsored  by  the 
Air  Force  Cambridge  Research  Laboratories  (AFSC) , United  States  Air 
Force,  under  contract  No.:  F44620-75-C-0029 . The  United  States  Govern- 
ment is  authorized  to  reproduce  and  distribute  reprints  for  govern- 
mental purposes  notwithstanding  any  copyright  hereon. 

The  experiments  referred  to  were  carried  out  by  A/S 
NORCONSULT,  Oslo  and  sponsored  by  the  Bundesministerium  der  Ver- 
teidigung.  West  Germany. 


'T 


d 


- 36  - 


REFERENCES 


[1]  Leif  N.  Persen,  "Rock  Dynamics  and  Geophysical  Exploration", 

Elsevier  Publ.  Co.,  Amsterdam  Oxford  New  York 
(1975),  p .53-67 . 

[2]  R.  Lenschow  and , "Carrying  Capacity  of  the  Intersection  Between 

Dome  and  Cylinder  Wall  of  a reinforced  Con- 
' crete  Structure",  A tentative  report  in  two 

sections,  presented  at  the  BOSS  '76  Conference 
("Behaviour  of  Offshore  Structures")  held  at 
The  Norwegian  Institute  of  Technology  (NTH) 
August  1976. 


APPENDIX 


Listing  of  the  static  quantities  in  Fig. 30: 

N = normalized  axial  force  = N/p^R 

Q = normalized  shearing  force  = Q/pQR 

M = normalized  bending  moment  = M/p^R2 

H = normalized  horizontal  force  = H/PQR 

V = normalized  vertical  force  = V/p^R 

Load : 

d-d  faosQ  - eosB  ) 
r ro  o 
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Q 

M 
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0.01 06927 
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0.0194581 
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0 

N 

Q 

M 
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